Determination of polarized PDFs from a QCD analysis of inclusive and semi-inclusive 

Deep Inelastic Scattering data 



O 

O 

Q 



a,; 

(N 
> 

in 
o 

d 



X 



Elliot LeadeiB 

Imperial College, Prince Consort Road, London SW7 2BW, England. 

Aleksander V. SidorcAQ 
Bogoliubov Theoretical Laboratory Joint Institute for Nuclear Research 141980 Dubna, Russia. 

Dimiter B. Stameno\0 

Institute for Nuclear Research and Nuclear Energy 
Bulgarian Academy of Sciences Blvd. Tsarigradsko Chaussee 72, Sofia 1784, Bulgaria 

Abstract 

A new combined next to leading order QCD analysis of the polarized inclusive and semi-inclusive 
deep inelastic lepton-hadron scattering (DIS) data is presented. In contrast to previous combined 
analyses, the 1/Q 2 terms (kinematic - target mass corrections, and dynamic - higher twist correc- 
tions) in the expression for the nucleon spin structure function gi are taken into account. The new 
COMPASS data are included in the analysis. The impact of the semi-inclusive data on the polarized 
parton densities (PDFs) and on the higher twist corrections is discussed. The new results for the 
PDFs are compared to our (Leader, Sidorov, Stamenov) LSS'06 PDFs, obtained from the fit to 
the inclusive DIS data alone, and to those obtained from the de Florian, Sassot, Stratmann, and 
Vogelsang global analysis. 

PACS numbers: 13.60.Hb, 12.38.-t, 14.20.Dh 



I. INTRODUCTION 

Experiments on polarized inclusive deep inelastic 
lepton-hadron scattering (DIS), reactions of the type 
I + p — > I' + X with both polarized lepton and hadron, 
because of the nonexistence of neutrino data, can only, in 
principle yield information on the sum of quark and anti- 
quark parton densities i.e. information on the polarized 
densities Au + Au, Ad + Ad, As + As and AG. 

Information about the antiquark densities Am, Ad and 
the separate As and As strange densities thus has to be 
extracted from other reactions, notably polarized semi- 
inclusive lepton-hadron reactions (SIDIS) l+p — s- l' + h + 
X, where h is a detected hadron in the final state, or from 
semi-inclusive hadronic scattering (SIHS) like p + p — > 
h + X, involving polarized protons, and only possible at 
the RHIC collider at Brookhaven National Laboratory. 

QCD analyses of polarized DIS data, at next to lead- 
ing order accuracy (NLO), have been carried out for some 
decades (for more recent analyses see [H, 0), but it was 
only in 2008 that de Florian, Sassot , Strattmann and 
Vogelsang (DSSV), in a ground-breaking paper Q, car- 
ried out a combined analysis of polarized DIS, SIDIS and 
SIHS, at NLO accuracy. 

The technical problems involved in going from an anal- 
ysis of DIS to such a combined analysis are formidable. 
In this paper we extend our study of polarized DIS to a 
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joint analysis of the world data on DIS and SIDIS reac- 
tions. 

In contrast to the situation in unpolarized DIS, a large 
portion of the most accurate data on polarized DIS lie 
in a kinematical region where target mass corrections 
(TMC) of order M 2 /Q 2 ( whose form is exactly known), 
and dynamical higher twist (HT) corrections of order 
Aq CD /Q 2 are important. We have thus included such 
terms in our description of the DIS data. However, for 
the SIDIS data, we do not know the analogous results 
at present, so do not include such terms. As it happens 
almost all the SIDIS data we utilize are in kinematic re- 
gions where such corrections should not be important. 

Despite the fact that it has been emphasized in the lit- 
erature for more than a decade that DIS data can only, in 
principle, yield information on the sum of quark and an- 
tiquark densities, some analyses of purely inclusive DIS 
continue to show results for valence densities, under what 
are termed assumptions about the sea-quark densities 
Au(x) and Ad(x). It is important to realize that these 
are not really physical assumptions, but merely conven- 
tions. In contrast, it is important to note that although 
we tend to think of the strange quark density as a sca- 
quark density, As(ir) + As(x) is fully determined by the 
purely inclusive DIS data. This is particularly important 
because of the apparent incompatibility of the polarized 
strange quark density obtained from DIS and from SIDIS 
data, as will be discussed in detail later. 

In this paper we present the results of our NLO QCD 
analysis of polarized inclusive and semi-inclusive DIS 
data. Our analysis differs from DSSV in the following 
respects: 

(i) We have included new data from the COMPASS 
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group at CERN, which were not available in 2008. 

(ii) We are more careful in handling the kinemat- 
ics and include target mass corrections and higher twist 
terms in the DIS sector of our analysis. 

(hi) Our parametrization of the parton densities is sim- 
ilar to that of DSSV, but differs in some important as- 
pects, as will be explained in detail in Sec. III. 



II. QCD FRAMEWORK FOR INCLUSIVE AND 
SEMI-INCLUSIVE POLARIZED DIS 



A. Inclusive DIS 

One of the features of polarized DIS is that more than 
half of the present data are at moderate Q 2 and W 2 
(Q 2 ~ 1 - 5 GeV 2 , 4 GeV 2 < W 2 < 10 GeV 2 ), or in the 
so-called preasymptotic region. This is especially the case 
for the very precise experiments performed at the Jeffer- 
son Laboratory. So, in contrast to the unpolarized case 
this region cannot be excluded from the analysis. As was 
shown in [4|, to confront correctly the QCD predictions 
to the experimental data including the preasymptotic re- 
gion, the nonperturbative higher twist (powers in 1/Q 2 ) 
corrections to the nucleon spin structure functions have 
to be taken into account too. 

In QCD the spin structure function g\ has the following 
form for Q 2 >> A 2 (the nucleon target label N is not 
shown): 



3i (x, Q 2 ) = 3i(a;,<3 2 ) L T + gi(x, Q 2 )ht 



(1) 



where "LT" denotes the leading twist (r = 2) contri- 
bution to 31, while "HT" denotes the contribution to 
3i arising from QCD operators of higher twist, namely 
t > 3: 

9i(x,Q 2 )lt = 3i(^Q 2 ) P QCD + ^ TMC (z,Q 2 )/Q 2 

+0(M 4 /Q 4 ) , (2) 

where gi(x, Q 2 ) p qcd is the well-known (logarithmic in 
Q 2 ) NLO perturbative QCD contribution 
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and h™ c (x,Q 2 ) are the exactly calculable kinematic 
target mass corrections [5J, which, being purely kine- 
matic, effectively belong to the LT term. In Eq. 
©, Aq(x,Q 2 ), Aq(x,Q 2 ) and AG(x,Q 2 ) are quark, an- 
tiquark and gluon polarized densities in the proton, 
which evolve in Q 2 according to the spin-dependent 
NLO DGLAP equations. 5C{x) q ^ G are the NLO spin- 
dependent Wilson coefficient functions and the symbol 
® denotes the usual convolution in Bjorken x space, nf 
is the number of active flavors (nj = 3 in our analysis). 



In addition to the LT contribution, the dynamical higher 
twist effects 

3i (x, Q 2 )ht = h(x, Q 2 )/Q 2 + 0(A 4 /Q 4 ) , (4) 

must be taken into account at low Q 2 . The latter are non- 
perturbative effects and cannot be calculated in a model 
independent way. That is why we prefer to extract them 
directly from the experimental data. The method used 
to extract simultaneously the polarized parton densities 
and higher twist corrections to 31 from data on gx/Fi 
and Ai(fa gx/Fi), is described in [J]. Accordingly the 
gi /F\ data have been fitted using the experimental data 
for the unpolarized structure function F\(x, Q 2 )'. 



gi(x,Q 2 ) 



Fi(x,Q 2 ) 



LT 
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(5) 



Jexp 



As usual, F\ is replaced by its expression in terms 
of the usually extracted from unpolarized DIS experi- 
ments i*2 and R. As in our previous analyses, the phe- 
nomenological parametrizations of the experimental data 
for F 2 (x,Q 2 ) [y] and the ratio R(x, Q 2 ) of the longitu- 
dinal to transverse 7^ cross-sections Q are used. Note 
that such a procedure is equivalent to a fit to (gi) e xp, but 
it is more precise than the fit to the 31 data themselves 
actually presented by the experimental groups because 
here the 31 data are extracted in the same way for all of 
the data sets. Note also, that in our analysis the loga- 
rithmic Q 2 dependence of h(x,Q 2 ) in Eq. (O, which is 
not known in QCD, is neglected. Compared to the prin- 
cipal 1/Q 2 dependence it is expected to be small and the 
accuracy of the present data does not allow its determi- 
nation. Therefore, the extracted from the data values of 
h(x) correspond to the mean Q 2 for each z-bean. 



B. Semi-inclusive DIS 

As in the inclusive DIS case, the measured double spin 
asymmetries A^ N for the polarized semi-inclusive deep 
inelastic scattering, I + N —> I + h + X, where in addition 
to the scattered lepton, hadron h is also detected, can be 
presented by the ratio of the spin structure functions g^ N 
and 32W' anc ^ ^ ne unpolarized structure function F^ N , 

A ^ X > Z > Q) - F? N (x,z,Q 2 ) ' (6) 

where x is the Bjorken variable, z — (ph-PN)/{pN-q) is 
the fractional energy of the hadrons in the the center-of- 
mass system (c.m.s) frame of the nucleon and the virtual 
photon, and q is the usual notation for the photon four- 
momentum (— q 2 = Q 2 ). In ^ the index N is used for 
the different targets and in what follows it will be sup- 
pressed. Note also that in ([6]) the contribution of the 
spin structure function 3/? to the asymmetry A\ can be 
neglected by two reasons. First, although g\ is not mea- 
sured yet, it is expected to be small as in the inclusive 
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DIS case. Second, the term is multiplied by a factor 
7 2 = 4M 2 x 2 /Q 2 which in the kinematic region of the 
present SIDIS experiments is negligible. 

For the time being it is not known how to account for 
the HT and TMC corrections in SIDIS processes. For- 
tunately, they should be less important due to the kine- 
matic region and the accuracy of the present SIDIS data. 
So, in our QCD analysis we will use the approximate 
equation: 



Ah / „ n 2N _ 9in( x > z >Q 2 )nlo 
F{ l N (x,z,Q 2 ) NLO 



(7) 



In NLO QCD the structure functions g\ and F± have 
the following forms: 



2g1(x,z,Q 2 ) 
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C^(a;, z) needed in these calculations can been found in 



For the unpolarized parton densities we use the NLO 
MRST'02 PDFs pi], and for the fragmentation func- 
tions, the NLO DSS set [Tljj for pions, kaons and unin- 
dcntificd hadrons. The main reason to use the MRST'02 
set for the unpolarized PDFs is that the DSS fragmenta- 
tion functions were extracted from the data using in the 
SIDIS sector the MRST'02 unpolarized PDFs and the 
corresponding a s (M§) value. 

Compared with our previous fits to the inclusive DIS 
data only (for example, see (2J) we use now a more gen- 
eral parametrization for the input NLO polarized parton 
densities at Q% = 1 GeV 2 . It has the form for (Ait + Aw) 
and (Ad + Ad) 



x(Au + Au)(x,Ql) 



A u+a x a ^(l 



(1 + e u+u \/x + ju+ax), 

x(Ad + Ad)(x,Ql) = A d+d x a ^(l-xf^ 
( l + £d+dVx + ^d+dx), 



(10) 



and 



xAu(x, Qq) 
xAd(x,Ql) 
xAs(x, Ql) 
xAG(x,Q 2 ) 



A a x a ^(l-x)^(l + j u x), 
A s x<**(l-xf*(l + ls x), 
A s x as (l-x)^(l+ Ts x), 
A G x aG (l~xf G (l + lG x), (11) 



for the polarized sea quarks Aq and the gluon parton 
densities. Since the accuracy of the present SIDIS data 
is not enough to distinguish As from As, we assume the 
relation As(x) = As(x). 

As usual, the set of free parameters {a^} in (fTU|) and 
(|11[) is reduced by the well-known sum rules 



a 3 = gA = F + D = 1.269 ± 0.003, 



(12) 



In Eqs. © and © AC\]\x,z) and C§\x,z) arc the 



= 3F-D = 0.585 ± 0.025, 



(13) 



NLO partonic coefficient functions in the MS scheme col- 
lected in Q. Djj, Dg are the fragmentation functions 
(FFs) for quarks, antiquarks and gluons, and nj is the 
number of active flavors (n/ = 3 in our present analysis). 



C. Method of analysis 

In our previous analyses of the inclusive DIS data the 
inverse Mellin transformation method has been used to 
calculate the spin structure function gi(x,Q 2 ) from its 
moments. The double Mellin transform technique was 
developed by Stratmann and Vogelsang and first applied 
in the NLO QCD analysis of SIDIS data 0. We have 
used it to calculate the structure functions gi(x,Q 2 ) 
and Fi(x,Q 2 ) from their moments. The expressions for 
the moments of the coefficient functions AC|.- (x, z) and 



where a 3 and as are nonsinglet combinations of the first 
moments of the polarized parton densities corresponding 
to 3 rd and 8 th components of the axial vector Cabibbo 
current, 



(Ad + Ad)(Q 2 ), 
(Ad + Ad)(Q 2 ) 



(14) 



(15) 



a 3 = (Au + Au)(Q 2 ) ~ 
a 8 = (Au + Au){Q 2 ) + 
-2(As + As)(Q 2 ) 



The constants gA in Eq. (fT2"|) and as in Eq. (fl"3"|) are 
taken from [l^] and QT| , respectively. The sum rule (fT2")) 
reflects isospin SU(2) symmetry, whereas (fTBl is a conse- 
quence of the SU (3) / flavor symmetry treatment of the 
hyperon /3-decays. So, using the constraints (|T2"j) and (fl"3"|) 
the parameters A u+a and A d+d in (1101) can be determined 
as functions of the rest of the parameters connected with 
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(Au + Au), (Ad + Ad) and As. In addition, we assume 
that the parameters a u+il and a d+( [ which control the 
small x behavior of (Au + Au) and (Ad + Ad) are equal 
to those of Au and Ad, respectively, which is reasonable 
as sea quarks likely dominate in this region. 

The large x behavior of the polarized sea quarks and 
gluon densities is mainly determined from the positivity 
constraints 

\Afi(x,Q 2 )\ < fi(x,Ql), \Afi(x,Ql)\ < Ux,Ql). 

(16) 

The constraints (fT6|) are the consequence of a proba- 
bilistic interpretation of the parton densities in the naive 
parton model, which is still valid in LO QCD. Beyond LO 
the parton densities are not physical quantities and the 
positivity constraints on the polarized parton densities 
are more complicated. They follow from the positivity 
condition for the polarized lepton-hadron cross-sections 
A(Ti in terms of the unpolarized ones (|A<7i| < Oi) and in- 
clude also the Wilson coefficient functions. It was shown 
[i~4j . however, that for all practical purposes it is enough, 
at the present stage, to consider LO positivity bounds 
for LO as well as for NLO parton densities, since NLO 
corrections are only relevant at the level of accuracy of a 
few percent. 

For the unpolarized NLO parton densities on the RHS 
of (Unj), we are using the MRST'02 parton densities [Icj . 
In order to guarantee fulfilling of the positivity condition, 
we assume the following equation for the parameters fa 
which control the large x behavior of the polarized sea 
quarks and gluons: 

fa = f3 G = p sea{M RST02) = 7.276. (17) 

The rest of parameters {a^} = {Ai, Oi, fa, ei,ji}, as 
well as the unknown higher twist corrections h N (x) to 
in (J5j) have been determined from a simultaneous fit to 
the DIS and SIDIS data. For the determination of HT the 
measured x region has been split into 5 bins and to any 
x-bin two parameters h\ and h\ have been attached 
For a deuteron target the relation h\ d) = 0.925(/i 4 (p) + 
h\ )/2 has been used. So, to the set of parameters {ai} 
connected with the input polarized PDFs (|T0l [TTj) . 10 

parameters for the HT corrections, hff* and h\ (i = 
1, 2, .., 5), have been added. 

In the polarized DIS and SIDIS processes the Q 2 range 
and the accuracy of the data are much smaller than that 
in the unpolarized case. That is why, in all calculations 
we have used a fixed value of the NLO QCD parame- 
ter A^jg(n/ = 4) = 311 MeV, which corresponds to 
a s (M 2 ) = 0.1197, as obtained by the MRST NLO QCD 
analysis [l(| of the world unpolarized data. The value of 
Aj^g above is slightly changed from that of MRST'02 be- 
cause the scale dependence of the strong running coupling 
a s (Q 2 ) is calculated using the so-called "Denominator" 



representation [15[ 

a s (Q 2 ) = 1 

^ ~ /3oMQ7Ay + ^IniHQ 2 /^) + %} ' 

which is a more precise iterative solution of its renor- 
malization group equation at NLO accuracy. In (|18[) 
/3 = 11 - 2n//3, fa = 102 - 38n//3 and A m (n f = 
3,4,5) = 366, 311, 224 MeV. The number of active fla- 
vors nf in a s (Q 2 ) was fixed by the number of quarks with 
ra\ < Q 2 taking m c = 1.43 GeV and m b = 4.3 GeV. 

The advantage of the analytic expression (fT5|) for a s is 
that: first, it reproduces with a very good accuracy the 
numerical solution of the renormalization group equation 
needed at small Q 2 , down to Q 2 = 1 GeV 2 , and second, 
for Q 2 > 4 GeV 2 it practically coincides with the be- 
havior of a s corresponding to its usual l/ln(Q 2 /A 2 -^-)- 
expansion at NLO (l^ |. 



III. RESULTS OF ANALYSIS 

The numerical results of our global NLO QCD fit to 
the world inclusive [l6l - [27| and semi-inclusive [2|| [28l - [3fl | 
DIS data are presented in Tables I and II. The data used 
(841 experimental points for DIS and 202 experimental 
points for SIDIS) cover the following kinematic regions: 
{0.005 < x < 0.75, 1 < Q 2 < 62 GeV 2 } for DIS and 
{0.005 < x < 0.48, 1 < Q 2 < 60 GeV 2 } for SIDIS 
processes. 

In our analysis the minimization of the \ 2 function is 
performed using the program MINUIT at CERN [3ll ] . 
The experimental errors are given by statistical and 
point-to-point systematic errors added in quadrature. In 
the minimum of \ 2 & n accurate (a positive definite) error 
matrix is obtained and the error bands of the polarized 
PDFs were calculated using the standard Hessian method 
with A\ 2 — 1. We understand that A% 2 = 1 could be 
an underestimation of the uncertainties of the polarized 
PDFs. If one wishes to use the choice A% 2 > 1, one has 
simply to scale our uncertainties of the polarized PDFs 
by (Ax 2 ) 1 / 2 - However, one has to keep in mind the fol- 
lowing points: 

(i) The systematic errors are partly correlated which 
would lead to an overestimation of the errors when added 
in quadrature with the statistical ones, which compen- 
sates a part of any underestimation arising from using 
A X 2 = 1. 

(ii) In the analyses of the groups which present un- 
certainties of polarized PDFs corresponding to A% 2 > 1 
(for an example, A^ 2 = UP ~ N, where N is the num- 
ber of the free parameters), only the statistical errors are 
usually taken into account. 

(iii) The different status of the PDFs and higher twist 
parameters - practically they are not correlated. 

(iv) The experimental errors in the polarized case are 
much larger then those in the unpolarized one. 
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TABLE I: Data used in our global NLO QCD analysis, the 
individual x 2 for each set and the total x 2 of the fit 



Experiment 


Process 


Ndata 


x 2 


EMC [16] 


DIS(p) 


10 


4.2 


SMC [17] 


DIS(p) 


12 


5.5 


SMC [17] 


DIS(d) 


12 


18.0 


COMPASS [18] 


DIS(p) 


15 


12.0 


COMPASS [19] 


DIS(d) 


15 


8.4 


SLAC/E142 [20] 


DIS(n) 


8 


5.8 


SLAC/E143 [21] 


DIS(p) 


28 


17.8 


SLAC/E143 [21] 


DIS(d) 


28 


39.9 


SLAC/E154 [22] 


DIS(n) 


11 


2.6 


SLAC/E155 [23] 


DIS(p) 


24 


25.5 


SLAC/E155 [24] 


DIS(d) 


24 


16.5 


HERMES [25] 


DIS(p) 


9 


5.4 


HERMES [25J 


DIS(d) 


9 


6.8 


JLab-Hall A [26] 


DIS(n) 


3 


0.3 


CLAS [27] 


DIS(p) 


151 


119.9 


CLAS [27] 


DIS(d) 


482 


427.9 


SMC [28] 


SIDIS(p,/i+) 


12 


18.1 


SMC [28] 


SIDIS(p,/i-) 


12 


11.2 


SMC [28J 


SIDIS(d,/i+) 


12 


9.4 


SMC [28J 


SIDIS(d,/i") 


12 


13.8 


HERMES [25] 


SIDIS(p,/i+) 


9 


5.9 


HERMES [25J 


SIDIS(p,/i-) 


9 


5.3 


HERMES [25] 


SIDIS(d,/i+) 


9 


10.5 


HERMES [25] 


SIDIS(d,/i") 


9 


4.8 


HERMES [25J 


SIDIS(p,7T + ) 


9 


9.9 


HERMES [25] 


SIDIS(p,7T ) 


9 


5.1 


HERMES [25] 


SIDIS(d,7T+) 


9 


8.6 


HERMES [25] 


SIDIS(d,7r ) 


9 


19.8 


HERMES [25] 


SIDIS(d,is: + ) 


9 


6.7 


HERMES [25] 


SIDIS(d,if~) 


9 


5.6 


COMPASS [29J 


SIDIS(d,/i+) 


12 


7.6 


COMPASS [29J 


SIDIS(d,/i") 


12 


10.9 


COMPASS [30J 


SIDIS(d,7T+) 


10 


2.6 


COMPASS [30J 


SIDIS(d,7r ) 


10 


4.5 


COMPASS [30J 


SIDIS(d,if+) 


10 


7.8 


COMPASS [30] 


SIDIS(d,Js-) 


10 


13.7 


TOTAL: 




1043 


898.6 




FIG. 1: Comparison of our NLO QCD results for the SIDIS 
asymmetries with the data at measured x and Q 2 . 



parameters e d+ j and 7j in Eqs. (|T0| and (fTTj) , respec- 
tively, so they were eliminated i.e. put equal to zero. 
Note that the central value of 7^ obtained from the fits 
was always positive. So that its elimination does not 
change the negative behavior of xAd(x) for any x in the 
measured region. 



In Table I the data sets, both for inclusive and semi- 
inclusive DIS, used in our analysis are listed and the cor- 
responding values of x 2 obtained from the best fit to 
the data are presented. As seen from Table I, a good 
description of the data is achieved for both the inclu- 
sive (XAr r p=0-85) and semi-inclusive (x% r p=0-90) pro- 
cesses (NrP is the number of corresponding experimental 
points). The total value of x 2 df ^ s 0.88. The quality of 
the fit to the data is demonstrated in Fig. 1 for some of 
the SIDIS asymmetries obtained by the HERMES and 
COMPASS Collaborations. 

The values of the parameters attached to the input po- 
larized PDFs obtained from the best fit to the data are 
presented in Table II. The errors correspond to A% 2 = 1 . 
Note also that only the experimental errors (statistical 
and systematic) are taken into account in their calcu- 
lation. It was impossible to determine from the fit the 



A. The role of semi-inclusive DIS data in 
determining the polarized sea-quark densities 

Let us discuss the impact of semi-inclusive DIS data 
on the polarized PDFs. Because of SIDIS data a flavor 
decomposition of the polarized sea is achieved and the 
light antiquark polarized densities Au(ir) and Ad(x) are 
determined without any additional assumptions. While 
Ad(x) is negative for any x in the measured x region, 
Au(x) is a positive, passes zero around x = 0.2 and 
becomes negative for large x. Sign-changing solutions 
are also found for the polarized strange sea As(x) and 
gluon AG(x) densities. The sign-changing behavior for 
AG(x) is not surprising since it was already found from 
the analysis of the inclusive DIS data alone 0. On the 
other hand, on the basis of results from all published 
analyses of inclusive DIS, we consider the sign-changing 
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TABLE II: The parameters of the NLO(MS) input polarized PDFs at Q 2 = 1 GeV 2 obtained from the best fit to the data. 
The parameters marked by (*) are fixed. 



flavor 


A 


a 


P 


e 


7 


u + u 


1.097* 


0.782 ± 0.165 


3.335 ± 0.154 


-1.779 ± 0.896 


10.20 ± 5.61 


d + d 


-0.394* 


0.547 ± 0.118 


4.056 ± 0.543 





6.758 ± 5.366 


u 


0.334 ± 0.174 


0.782 ± 0.165 


7.267* 





-5.136 ± 2.414 


(1 


-0.133 ± 0.075 


0.547 ± 0.118 


7.267* 








s 


-0.00352 ± 0.00194 


0.0458 ± 0.0206 


7.267* 





-39.02 ± 29.62 


G 


-68.23 ± 65.79 


1.975 ± 0.459 


7.267* 





-3.536 ± 1.089 



solution for As(x) quite puzzling. The central values 
of the sea-quark and gluon polarized densities together 
with their error bands are presented and compared to 
those of DSSV (dashed curves) in Fig. 2. The main dif- 
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Q 2 = 2.5 GeV 2 


LSS'10 
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0.1 1 
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2: Our NLO sea quarks and gluon polarized PDFs at 
2.5 GeV 2 in the MS scheme. For comparison the DSSV 
PDFs 0] are also presented. 



ference between the Leader-Sidorov-Stamenov (LSS) and 
DSSV sets is in the strange sea-quark density As(x). Al- 
though the first moments are almost equal (-0.054 and 
-0.055 at Q 2 = 1 GeV 2 for LSS and DSSV, respectively), 
our As(x) is less negative for x < 0.03 and less positive 
for x > 0.03. Note that DSSV have used an additional 
constraint a s = aj for the parameters ag and aj which 
means a similar small x behavior for the sea-quark den- 
sities As(x) andAd(x). We do not think this assumption 
is reasonable. The central values of our gluon density 
and its first moment are rather different from those of 
DSSV, however they coincide within the errors which are 
still large in the measured x region. 

In Fig. 2 our LSS'06 PDFs (dotted curves) 0] obtained 
from the NLO QCD analysis of the world inclusive DIS 
data are presented too. While the light antiquark polar- 
ized densities Au(x) and Ad(x) cannot be, in principle, 
determined from polarized inclusive DIS data, the sum 
(As + As)(x, Q 2 ) is well determined and all the NLO 



QCD analyses yield for this sum a negative value for any x 
in the measured region (for example, see Refs. [EI 11). 
In these analyses, however, a term like (1 + jx), which 
would permit a sign-change, was not included in the input 
parametrization of (As + As)(x, Qq)/2 (33[. We therefore 
re-analyzed the world polarized inclusive DIS data using 
such a term in the input strange sea-quark density 

(As + As)(x, Ql)/2 = Ax a (l - xf(l + jx). (19) 

Our preliminary results confirm the previous ones, 
namely, that (As + As)(x, Q 2 )/2 is negative in the mea- 
sured (x, Q 2 ) region. So, the behavior of the polarized 
strange quark density remains controversial. Note that 
in the presence of SIDIS data As and As can, in prin- 
ciple, be separately determined, as was done recently 
by the COMPASS Collaboration, where it was shown 
[34| that there is no significant difference between As(x) 
and As(x) in the x-range covered by their inclusive and 
semi-inclusive DIS data. This result was obtained in the 
LO QCD approximation, with the additional assumption 
that the SIDIS asymmetries are Q 2 -independent. We 
checked the latter assumption using in the calculations 
of the asymmetries our NLO PDFs, and found it not 
quite correct. Also, the errors of the extracted values of 
the difference x(As(x) — As(x)) are rather too large to 
allow us to conclude that the assumption As(x) = As(x) 
used in our analysis and that of the DSSV is correct. So, 
if it is not correct, it might possibly be the cause that 
(As + As)(x, Q 2 )/2 densities obtained from the analyses 
of inclusive DIS data and combined inclusive and semi- 
inclusive DIS data sets, respectively, are in contradiction. 
However, at first glance, it looks as if the difference be- 
tween As and As would have to be quite significant and 
might contradict the COMPASS results. Perhaps a more 
important issue is the sensitivity of the results to the 
form of the fragmentation functions. An analysis by the 
COMPASS group [3(3] demonstrated that the determina- 
tion of As(x) strongly depends on the set of the fragmen- 
tation functions used in the analysis and that the DSS 
FFs are crucially responsible for the unexpected behav- 
ior of As(x) obtained from the combined analysis. The 
study of the sensitivity of As(x) to different sets of FFs 
used in the analysis is one of the key points we plan to 
investigate in the future. 

In Fig. 3 we present our results for the polarized Au(x) 
and Ad(x) densities at Q 2 = 2.5 GeV 2 , which are con- 
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FIG. 3: Our results for Aw, Ad, (Aw + Au) and (Ad + Ad) 
polarized parton densities at NLO approximation. DSSV 
as well as LSS'06 2] results for the corresponding densities 
are presented too. 



TABLE III: The values of higher twist corrections extracted 
from the data in a model independent way. 





h p (xi) [GeV 2 ] 




h n (xi) [GeV 2 ] 


0.028 


-0.048 ± 0.037 


0.028 


0.093 ± 0.051 


0.100 


-0.084 ± 0.016 


0.100 


0.041 ± 0.036 


0.200 


-0.053 ± 0.012 


0.200 


0.000 ± 0.027 


0.350 


-0.045 ± 0.011 


0.325 


0.022 ± 0.018 


0.600 


-0.012 ± 0.014 


0.500 


0.017 ± 0.015 



sistent with those obtained by DSSV (dashed curves). 
As expected, the SIDIS data do not influence essentially 
the sums (Au(x) + Au(x)) and (Ad(x) + Ad(x)) already 
well determined from the analysis of the inclusive DIS 
data. This fact is illustrated in Fig. 3 where our re- 
sults from the combined analysis are compared with our 
LSS'06 PDFs. 



B. High twist effects 

As mentioned in Sec. II A, in contrast to other com- 
bined analyses of the inclusive and semi-inclusive DIS 
data, we take into account the target mass and higher 
twist corrections in the DIS sector. The values of the HT 
corrections to gi extracted from the data in this analysis 
are presented in Table III and illustrated in Fig. 4. 

Compared to the HT(LSS'06) corrections obtained in 
our analysis of the inclusive DIS data alone jlj the values 
of the HT corrections for the proton target are practically 
not changed, while the central values of HT corrections 
for the neutron target are smaller in the region x < 0.2, 
but in agreement with HT^ 11 ) (LSS'06) within the errors, 
excepting the x region around x = 0.1. We consider 
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FIG. 4: Comparison between HT values corresponding to the 
fits of inclusive DIS @] and combined inclusive and SIDIS 
data set (this analysis). 



the tendency of the HT^ 11 ) corrections to be smaller in 
the region x < 0.2 to be a result of the new behavior 
of As(x), i.e. positive for x > 0.03. The positive con- 
tribution in <?" from As(x) should be compensated by a 
less positive HT^ n ) contribution in this region. Since the 
biggest difference between the values of As(:t)(dis+sidis) 
and As(x)dis is in the region x ~ 0.1 (see Fig. 2), this 
effect is biggest in this x region. The impact of As(x) 
on HT corrections is visible mainly for the neutron tar- 
get because the contribution of As(x) in g\ is relatively 
larger than that in g p . 

Let us briefly discuss the values of the first moments 
of the higher twist corrections to the proton and neu- 
tron structure function g\. Using the values for h N (x) 
from Table III we obtain for their first moments in the 
experimental region: 



/■0.75 

h N = / 

JO. 0045 



iV 



(x)dx, (N — p, n) 



(20) 



h p = (-0.028 ± 0.005) GeV 2 for the proton and h n = 
(0.018 ±0.008) GeV 2 for_the neutron target. As a result, 
for the nonsinglet (h p — h n ) and the singlet (h p + h n ) we 
obtain (-0.046±0.009) GeV 2 and (-0.011±0.009) GeV 2 , 
respectively. The statistical and systematic errors are 
added in quadrature. Note that in our notation h — 
f Q h(x)dx = 4M 2 (d2+/2)/9, where di and /2 are the well 
known quantities, connected with the matrix elements of 
twist 3 and twist 4 operators, respectively [35j . and M is 
the mass of the nucleon. 

Our value for the nonsinglet (h p — h n ) is well consistent 
within the errors with those extracted directly from the 
recent analyses of the first moment of the nonsinglet spin 
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structure function 



(v- 



,Q 2 ) 



Note that our 



value for the nonsinglet (h p — h n ) is also in agreement 

as well as with 
The values ob- 



38 
40 



with the QCD sum rule estimates 
the instanton model predictions [39 
tained for the nonsinglet (hP — h n ) and singlet (h p + h n ) 
quantities are in qualitative agreement with the relation 
\h p + « \hP - h n \ derived in the large N c limit in 
QCD H3. 

Our results on the higher twist effects are not in agree- 
ment with those obtained in (4l| , where the authors find 
no evidence for HT. 



C. The sign of the gluon polarization 

We have found that the combined NLO QCD analy- 
sis of the present polarized inclusive DIS and SIDIS data 
cannot rule out the solution with a positive gluon polar- 
ization. The values of \ 2 /DF corresponding to the fits 
with sign-changing xAG(x, Q 2 ) and positive xAG(x, Q 2 ) 
are practically the same: x 2 /DF (node xAG) — 0.883 
and x 2 / DF (xAG > 0) = 0.888, and the data cannot 
distinguish between these two solutions (see Fig. 5). 
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FIG. 5: Comparison between the positive and sign-changing 
gluon densities. The corresponding error bands are also 
shown. 

The corresponding sea-quark densities are shown in 
Fig. 6. As seen, the sea-quark densities obtained in the 
fits with positive and sign-changing xAG(x) are almost 
identical. Note that the extracted HT values correspond- 
ing to both fits are also effectively identical. As a result, 
one can conclude that including the SIDIS data in the 
QCD analysis does not help to constrain better the po- 
larized gluon density. 

In Fig. 7 the ratio AG(x)/G{x) calculated for both the 
sign-changing and positive solutions for AG(x) obtained 
in our NLO QCD analysis is compared with the directly 
measured values of AG/G obtained from a quasireal pho- 
toproduction of high p t hadron pairs |42h44| ■ and from 
the open charm production (45j measurements. For the 
unpolarized gluon density G{x) in the ratio above we 
have used that of the NLO MRST'02 0. The theoret- 
ical curves are given for /i 2 = 3 GeV 2 (high p t hadron 



FIG. 6: Comparison between the sea-quark densities corre- 
sponding to positive and changing in sign gluon densities. 



pairs) and fj 2 = 13 GeV 2 (open charm). As seen from 
Fig. 7, both solutions for the polarized gluon density are 
well consistent with the experimental values of AG/G. 
It should be noted, however, that in the extraction of 
AG/G by the experiments a LO QCD treatment has 
been used. A NLO extraction of the measured values 
is needed in order for this comparison to be quite cor- 
rect. In conclusion, the magnitude of the gluon density 
xAG(x) obtained from our combined NLO QCD analy- 
sis of inclusive and semi-inclusive DIS data and indepen- 
dently, from the photon-gluon fusion processes, is small 
in the region x ~ 0.08 — 0.2. 

When this analysis was finished, the COMPASS Col- 
laboration reported the first data on the asymmetries 
A^ + p~\ A 1 p <_) for charged pions and kaons produced 
on a proton target 34]. As seen in Fig. 8, our predic- 
tions for these asymmetries are in very good agreement 
with the data at measured x and Q 2 . 



D. The spin sum rule 

Let us finally discuss the present status of the proton 
spin sum rule. Using the values for AE(Q 2 ) and AG(Q 2 ) 
at Q 2 = 4 GeV 2 , the first moments of the quark singlet 
AS(x, Q 2 ) and gluon AG(x,Q 2 ) densities, obtained in 
our analysis (see Table IV) one finds for the spin of the 
proton: 
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TABLE IV: First moments of polarized PDFs at Q 2 = 4 GeV 2 . The corresponding DSSV values are also presented. 
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AE 


LSS10 (pos xAG) 


-0.063 ± 0.004 
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0.207 ± 0.034 
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-0.055 ± 0.006 
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FIG. 7: Comparison between the experimental data and 
NLO(MS) curves for the ratio AG{x)/G{x) at Q 2 = 3 GeV 2 
(top - high p t pairs) and Q 2 = 13 GeV 2 (bottom - open 
charm) corresponding to positive and sign-changing xAG. Er- 
ror bars represent the total (statistical and systematic) errors. 
The horizontal bar on each point shows the a;-range of the 
measurement. The NLO AAC (second listing of Ref. in []]]) 
and DSSV Q| curves on AG(x)/G(x) are also presented. 



Jz 



1 



A£(Q 2 )+AG(Q 2 ) + L Z (Q 2 ) 



= -0.21 ± 0.46 + L Z (Q 2 ) (node AG), 
= 0.42 ± 0.19 + L Z (Q 2 ) (pos AG). 



(21) 



In Eq. ([2~Tj) L Z (Q 2 ) is the sum of the angular orbital 
momenta of the quarks and gluons. Although the central 
values of the quark-gluon contribution in (|2ip are very 
different in the two cases, in view of the large uncertainty 
coming mainly from the gluons, one cannot yet come to a 
definite conclusion about the contribution of the orbital 
angular momentum to the total spin of the proton. 




FIG. 8: Our predictions for the COMPASS asymmetries for 
charged pions and kaons produced on a proton target. 



IV. SUMMARY 

A new combined NLO QCD analysis of the polar- 
ized inclusive and semi-inclusive DIS data is presented. 
In contrast to previous combined analyses, the 1/Q 2 
terms (kinematic - target mass corrections, and dynamic 
- higher twist corrections) to the nucleon spin structure 
function g\ are taken into account. The new results for 
the PDFs are compared to both the LSS'06 PDFs ob- 
tained from a fit to the inclusive DIS data alone, and to 
those obtained from the DSSV global analysis. The role 
of the semi-inclusive data in determining the polarized 
sea quarks is discussed. Because of SIDIS data Au(x, Q 2 ) 
and Ad(ir, Q 2 ), as well Au(x, Q 2 ) and Ad(x, Q 2 ) are de- 
termined without additional assumptions about the light 
sea quarks. The SIDIS data, analyzed under the assump- 
tion As(x,Q 2 ) — As(x 7 Q 2 ), imposes a sign-changing 
As(x, Q 2 ), as in the DSSV analysis, but our values are 
smaller in magnitude, less negative at x < 0.03 and less 
positive for x > 0.03. Note that As(x, Q 2 )sidis differs 
essentially from the negative ^(As + As)(x, Q 2 )dis ob- 
tained from all the QCD analyses of inclusive DIS data. 
As was mentioned above the behavior of As(x, Q 2 )sidis 
strongly depends on the fragmentation functions used in 
our analysis and that of the DSSV. A further detailed 
analysis of the sensitivity of As(x, Q 2 ) to the FFs is 
needed, and any model independent constraints on FFs 
would help. Another possible reason for this disagree- 
ment could be the assumption As(x, Q^) = As(x, Qq) 
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made in the global analyses. However, this would prob- 
ably require a significant difference between As and As. 
In any case, obtaining a final and unequivocal result for 
As(x) remains a challenge for further research on the 
internal spin structure of the nucleon. 

We have found also that the polarized gluon density 
is still ambiguous, and the present polarized DIS and 
SIDIS data cannot distinguish between the positive and 
sign-changing gluon densities AG(x). This ambiguity is 
the main reason that the quark-gluon contribution into 
the total spin of the proton is still not well determined. 

Finally, our combined NLO QCD analysis confirms our 
previous results on the higher twist corrections to the 



nucleon spin structure function , namely, that they are 
not negligible in the preasymptotic region and have to be 
accounted for in order to extract correctly the polarized 
PDFs. 
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